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Abstract 

A 3D nonlinear chaotic system, called the T system, is analyzed in this paper. Horseshoe 
chaos is investigated via the heteroclinic Shilnikov method constructing a heteroclinic connec- 
tions between the saddle equilibrium points of the system. Partially numerical computations 
are carried out to support the analytical results. 

1 Introduction 

One of the first examples of continuous dynamical system of dimension three, whose numerical 
simulations display the property of sensitivity to initial conditions, is the Lorenz system [j]]: 

x = a(y — x), y — cx — y — xz, z = —bz + xy, (1) 

where a, b, c are positive real parameters. Physical realization of the Lorenz system is the Rayleigh- 
Benard experiment. The system was derived from the hydrodynamical Navier-Stokes equations to 
create a dynamical model for meteorology. 

Trying to transform the Lorenz system from a stable state to a chaotic one (concept known 
as anticontrol of chaos or chaotification, 0] ), Chen in 3 introduces a 3D polynomial system 
(of type Lorenz) known as the Chen system. Vanecek and Celikovsky in |19j . characterize a 
generalized Lorenz system by a condition on its linear part matrix J = (ay): 0,120*21 > 0. Lii in 
[H] observes that while the classical Lorenz system belongs to this class, the Chen system satisfies 
a dual condition 012021 < 0, (systems that satisfy such a condition are called dual systems to the 
Lorenz system) and introduces a new system that bridges the gap between the Chen system and 
Lorenz system, satisfying 012021 = 0. A nonlinear system arising from a nuclear spin generator 
and compared with the Lorenz system is studied in |12| . Nonlinear dynamics is met in many areas 
from Ecology ^5] to Physics [7]. 

From the point of view of the potential applications, systems with sensitivity to the initial 
conditions can be used in secure communications pQ , |15| . Of the pioneering papers which proposed 
to use the chaotic systems in communications are the papers of Pecora and Carrol |10j . 
Consequently, an appropriate chaotic system can be chosen from a catalogue of chaotic systems to 
optimize some desirable factors, idea suggested in |15| . 

These some ideas led us to study a new 3D polynomial differential system given by: 

x = a{y — x), y — (c — a)x — axz, z = —bz + xy, (2) 
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with a, b, c real parameters and a/0. Call it the T system. Some results regarding the T system 
are already presented in J7| and ^H] . Compared with the Lii system introduced in jS] , the system 
T allows a larger possibility in choosing the parameters of the system and, therefore, it displays a 
more complex dynamics. 

The paper is organized as follows. In the second Section we give some details of the equilibria. 
In Section 3, we pay attention to the chaos in the system T and, using a method of type Shilnikov, 
we prove that the system displays "horseshoe" chaos. This implies that the system possesses a 
strange chaotic attractor. Not any strange dynamics is chaotical as it is pointed out, for example, 
in where is presented a scenario for a 2-dimensional dynamic system to possesses strange 
nonchaotic behavior. 

2 The equilibrium points of the system T 

Straightforward computations lead us to the following result: 

Proposition 2.1. If -(c — a) > 0, the system T possesses three equilibrium isolated points: 

and for b ^ 0, -(c — a) < it has only one isolated equilibrium point, O(0, 0, 0). 

As reported in |17|. we have: 
Theorem 2.1. For b ^ the following statements are true: 

a) If (a > 0, b > 0, c < a), then O(0, 0, 0) is asymptotically stable , 

b) lf(b<0) or (a < 0) or (a > 0, c > a), then O (0, 0, 0) is unstable. 

Theorem 2.2. If (a + b > 0, ab(c — a) > 0, b(2a 2 + bc — ac) > 0), the equilibrium points 
i?i.2(±Y / '^(c — a), ±^/^(c — a), are asymptotically stable. 

3 Chaos in the T system via the heteroclinic Shilnikov method 

Computing the Lyapunov exponents with the software Dynamics 9 for the parameter vector 
(a,b,c) = (2.1,0.6,30) and the initial conditions (0.1,-0.3,0.2), we get that Ai = 0.37 > 0, 
A2 = 0.00 and A3 = —3.07. So the system T displays chaotic characteristics. 

In the following, using the Shilnikov heteroclinic method, we show that the system T presents 
chaos of horseshoe type. We will employ the following result: 

Theorem 3.1. (\14\ If a 3-D given system x = F(x) has two equilibrium points Ex, E2, of 
type saddle-focus, i. e. the eigenvalues of the Jacobi matrix associated to the system in these points 
are 71 £l and ak ± ifik € C, fc = 1, 2, such that 

ai«2 > or 7172 > (3) 

and (the Shilnikov inequality) 

\ lk \ > |a fc |,fc = l,2, (4) 

and if the system has a heteroclinic orbit connecting the equilibrium points E\ and E%, then the 
Poincare map defined on a transversal section of the flow in a neighborhood of the heteroclinic orbit 
presents chaos of horseshoe type. 
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We will show that our system T given by 

x = a(y — x), y — (c — a)x — axz, z = —bz + xy, (5) 

fulfills the conditions from the above theorem. The equilibrium points of the system T are: 
Ei_2^t\J ^(c — a),±\J ^(c — a), for ^(c— a) > 0. The characteristic polynomial associated to 
the Jacobi matrix of the system T in these points is: 

/ (A) = A 3 + A 2 (a + b) + 6cA + 2ab(c - a) = (6) 

Denoting A = fi — (a + 6)/3, © leads to: 

f i 3 +PH + q = 0, (7) 

where 

p= be- -a 2 - -ab- -b 2 (8) 

^333 y ' 

and 

q = — a 3 a 2 b + -ab 2 H b 3 + -bca b 2 c. (9) 

y 27 9 9 27 3 3 w 

Denote A = (|) 2 + (|) 3 . Then, if A > 0, Eq. Q has a negative solution, a±, together with a 
pair of complex solutions, cti ± 10:3, where 



», = ,/-| + Va+«-|-Va, 



So, if A > 0, the three roots of Eq. JBJ are: 

. a+fe a+6 . a+6 . 

Ai = h ax, A 2 = h«2 + 1^3, A 3 = - ha 2 - ia 3 , (10) 

with Ai < 0. From the first two Eqs. of © one gets 
V = x + % 

z = - ^—^ = + and using the last Eq. of © we get 

ft {^tr) + b!t ir + a2x2 + axi -ab{c-a) = 0, 
or, equivalently 

x'x + (a + b) xx — xx — ax 2 + abxx + a 2 x A + ax 3 x — ab(c — a) x 2 = 0. (11) 
We observe that if x (t) is determined, then we can find y(t),z(t). Therefore, we have to find 
a function <p(t) such that x(t) — <p(t) to satisfy Eq. (|1 1|1 and <p(t) — > — \J ^(c — a) for t — > +00, 
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4>{t) — ► \J ^(c — a) for t — > — oo or, inversely 0(f) — > y^~(c — a) for < — * +00, — > — \J ^{c — a) 
for t — > —00. 

We can define, without loosing the generality, that the direction from £7j to _B 2 corresponds to 
f — > 00, and from E2 to £1 corresponds to t — > —00. 



00 



Consider <^(t) of the following form <j> (t) = —xq + a n e ™ Q * , with a real number and 



x o = y;(c-a). 

Identifying the coefficients of e" Qt in we get: 

x [a 3 + a 2 (a + b) + bca + 2ab(c - a)] a x = 0, (12) 

a 2 = t~, — v — (a 2 b + aba — ab(c — a) + 6a 2 XQ + iaaxV) a?, (13) 
/ (2a) x 



° 3 = ~ f(3a)x (4a 2 a;oai + 3aax af) + 

+ /(3a)a l^ a3 + a2fl + ^ct 2 & + 3a6a — 2a6c + 2a 2 b + \2a 2 x\ + 9aaxl] aia 2 , 
and for n > 4 

a « = 7(^( fl2 %9+ a M + 7(iH £ a 3 (j 3 - y 2 ) 0^+ 

i-\-j=n 

+ f(na)x„ E t" 2 (( a + fo ) J' 2 _ y°) + aba J - ab(c- a) + Qa 2 xl + iaajx 2 ] 

i-\-j—n 



(14) 



(15) 



where 

/ (na) = (na) 3 + (na) 2 (a + b) + bcna + 2ab(c — a), (16) 



ttijp, = -4xo 2J a i a ] a p + aia 3 a p a q , (17) 

i+j+p=n i+j'+p+g=n 

&iy P9 = — 3axo paidjCip + a pa,iaja p a q , (18) 

i+j+p=n i+j+p+q=n 

with i, j,p,<? > 1. 

In order to can be implemented on the computer, we put the sums (|17|l and (|18|l in the following 
forms: 

n — 2 / n—i—1 \ n — 3 / n—i—1 /n—i—j—1 

a vpq = _ 4xo^a.( ^ a.,a n _j_j J + ^ a, I ^ a 5 I ^ a p a n -i-j- v J , (19) 

i=i \ i=i y i=i \ j=i V P =i 

n — 2 / n—i—1 \ 



i=l \ 3=1 



Ti — 3 / ri — i— 1 / n—i—j — 1 
i=l \ ] = 1 V p=l 



(20) 
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Assume ai ^ 0. If not, we observe inductively that all coefficients a n ,n > 2, are zero. Hence, 
from JT^J, we get 

a 3 + a 2 (a + b) + bca + 2ab(c- a) = (21) 

that is a is the negative root of characteristic polynomial iJBJ. Because Eq. JHJl has a single negative 
solution for A > 0, we get that 

/ (not) = (naf + (naif (a + b) + bcna + 2ab(c - a) ^ 0,n > 1. (22) 

Consequently, the coefficients a n are completely determined by a,b,c,a and oi, and they are of 
the following form: 

a„ = g(n)a",n > 1 (23) 
where the terms g(n) are known functions. So the corresponding branch of the heteroclinic orbit 

oo 

for t > is determined. In a similar manner, consider the case t < and <j> (t) = xo + "ne""'", 

n=l 

with P real, and find (3 = a and b n = —a n , n > 0. Consequently, the heteroclinic orbit <p(t) is given 
by: 

y OO 

-xo + £ a n e nat for t>0 

71=1 

• (t) = { for t = (24) 

OO 

x - E a n e~ nat for t<0 

71=1 

Following a method presented in |21j one can show that <p(t) is uniformly convergent. So we have 
the following result that characterizes the chaos in the system T: 

Theorem 3.2. If A > and a\ + oi2 < — 2 (°+ fc ) ; then the system T has a heteroclinic orbit given 
by \24\ , which connects the equilibrium points Ei,E2, so the chaos is of horseshoe type. 

Let us partially illustrate the particular case a = 2.1, b = 0.6, c = 30. Then \± = a = —3.429, 
A 2 = 0.364-4.513i, A 3 = 0.364+4.513«, A = 911.69 > 0, x Q = 2.8234. Observe that the equilibrium 
points Ei, E 2 are saddle-focus. Imposing that <j>(t) to be at least continue, from </>(0_) = 0(0+) = 0, 
we find the equation 

— xo + ai + a 2 + ... + a n =0, n > 1. (25) 
Solving this equation, we observe that the first coefficient is a\ — 3.051 for any n > 10. 

Remark 3.1. Numerical series \2J$ which describes the heteroclinic orbit between E\ and E 2 , is 
rapidly convergent. For example, considering the first ten terms of the series, for t — 10 we get 

10 10 

<t>(t) = -xo + E a n e nat = -2.8234 7 for t = -10, <j>(t) = x - E a n e- nat = 2.8234 and for t = 0, 

n—l n—1 
10 

(j)(t) = —xo + E a n£ nat = 0.0000. So the equilibrium points E\, E 2 and O are found with an 



approximation of four exact decimals because xq = \J ~~~ a ) 



2.8234. 
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4 Conclusions 



In this paper we investigated a new chaotic system. For some parameter vectors, the system 
presents the property of dependence to initial conditions which is a necessary condition for a system 
to be chaotic. The system possesses three equilibrium points, the origin O(0, 0, 0) and another two 
points Ei_2- The chaos of horseshoe type was analyzed using the heteroclinic Shilnikov method, 
constructing a heteroclinic connection between the saddle- focus equilibrium points E\ and E^- 
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